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We show how to extract the scale-dependent mass anomalous dimension from the renormalization 
group invariant mode number of the massless Dirac operator. By combining simulations on multiple 
lattice volumes, and when possible several gauge couplings, we are able to measure the anomalous 
dimension across a wide range of energy scales. The method that we present is universal and can 
be applied to any lattice model of interest, including both conformal or chirally broken systems. 
We investigate SU(3) lattice gauge theories with Nf = 4, 8 and 12 light or massless fermions. The 
4-flavor model behaves as expected for a QCD-like system and demonstrates that systematic effects 
can be controlled in practical lattice calculations. Our 12-flavor results indicate the existence of 
an infrared fixed point, at which we predict the scheme-independent mass anomalous dimension 
7m — 0.32(3). For the 8-flavor model we observe a large anomalous dimension across a wide range 
of energy scales. Further investigation is required to determine whether Nf — 8 is chirally broken 
and walking, or if it possesses a strongly-coupled conformal fixed point. 

PACS numbers: 11. 15. Ha, 11.25.Hf, 12.60.Nz 



I. INTRODUCTION AND OVERVIEW 

The eigenvalues A of the fermion Dirac operator have 
been used extensively to gather information about the 
dynamics of lattice QCD. The most common application 
of Dirac eigenmodes has been the calculation of the chi- 
ral condensate E. In chirally broken systems, random 
matrix theory (RMT) predictions for the distributions of 
the lowest-energy modes in the e-regime allow the deter- 
mination of E at modest computational cost [1-4]. 

Another promising observable is the mode number 

i/(A) = V J p(u)duj, (1) 

where p(A) is the spectral density of the massless Dirac 
operator. Even though many more eigenmodes are 
needed to analyze the mode number ^(A), its renormal- 
ization group invar iance [5] makes it an ideal observable. 
Ref. [5] used a stochastic method to evaluate the mode 
number on a set of 2-flavor SU(3) p-regime configurations 
in an extended A range. The slope of v(\) predicts the 
spectral density, which is extrapolated to the chiral limit 
to determine E using the Banks-Casher relation [6]. 

Dirac eigenmodes have been considered in infrared- 
conformal gauge theories as well. In these systems the 
low-energy behavior is governed by a conformal infrared 
fixed point (IRFP) and there is no spontaneous chiral 
symmetry breaking. While one can search for IR con- 
formality by checking whether distributions of low-lying 
eigenmodes deviate from RMT predictions [7, 8] , this ap- 
proach is complicated by the lack of comparably rigorous 
theoretical predictions for eigenvalue distributions in IR- 
conformal systems. An alternative proposed by Ref. [9] 
is to investigate the finite-volume scaling of individual 
eigenmodes, which is related to the scheme-independent 
mass anomalous dimension 7^ at the IR fixed point. In 
Ref. [10] we extended this approach to describe the finite- 



volume scaling of several low-energy eigenmodes simul- 
taneously. In the time since the pilot study reported in 
Ref. [10], we have found the volume-scaling of low-energy 
eigenmodes to show fairly large systematic effects that 
are difficult to address in the absence of a more rigorous 
theoretical framework. 

In IR-conformal systems, the mode number appears 
to be a more robust quantity than individual eigen- 
modes. Based on the RG invariance of the mode number, 
Ref. [11] showed how the scaling of v(X) with A in the in- 
frared limit is related to the scheme-independent mass 
anomalous dimension at the IR fixed point. Ref. [11] 
also attempted to use this relation to extract 7^ for 
SU(2) gauge theory with two adjoint fermions, which 
is believed to be IR conformal [12-18] (see also the re- 
views [19, 20]). The initial results were plagued by large 
uncertainties, mainly because relatively few modes were 
calculated on small volumes (only the lowest 200 eigen- 
values on 16 3 x32 lattices) [11]. Ref. [21] used a stochastic 
method [5] to calculate the mode number v(\) across a 
much wider range of A, which allowed a stable fit, predict- 
ing 7^ = 0.371(20) for the same model. Similar studies 
have begun for SU(7V) gauge theories with two adjoint 
fermions, in the large- N limit [22]. 

In this paper we argue that the eigenvalues of the Dirac 
operator contain much more information than E or 7^, 
the targets of the calculations mentioned above. The 
scaling of the mode number v(\) depends on A, which 
defines an energy scale. We show how this behavior can 
be extracted from lattice simulations of both QCD-like 
and infrared-conformal systems, to determine the scale- 
dependent mass anomalous dimension 7 m (/i). In Ref. [23] 
we presented preliminary results of this approach, and 
demonstrated that combining data from several lattice 
volumes is an effective way to cover a wider range of 
energies. Here we improve this technique by also investi- 
gating multiple gauge couplings. For QCD-like systems 
we are even able to combine results at different gauge 
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FIG. 1. Our results for the mass anomalous dimension 7 m for SU(3) gauge theories with Nf = 4 (a and b, top), Nf = 8 (c, 
bottom left) and Nf — 12 (d, bottom right). In panel (b) we rescale the Nf — 4 results to be expressed in terms of a common 
lattice spacing (07.4 corresponding to /3f = 7.4). The displayed error bands account for error bars from a jackknife analysis as 
well as the fit ranges AA listed in Tables I through III. 



couplings, which extends the covered energy range and 
allows us to follow the evolution of the system from the 
perturbative UV limit to the onset of chiral symmetry 
breaking in the IR. 

The approach we propose is very general and can be 
used with any lattice model, offering a new way to in- 
vestigate the scale dependence of both IR-conformal and 
chirally broken systems. We review the energy depen- 
dence of the mass anomalous dimension and our method 
to extract it from lattice data in Section II, where we 
also discuss potential systematic effects. In particular, 
since we will present results from lattice calculations car- 
ried out with very light or massless fermions, we carefully 
consider the question of finite- volume effects. We show 
how combining different volumes allows us to ensure that 
our results correspond to the infinite-volume physics. 

We test our proposal in Sections III and IV for SU(3) 
lattice gauge theories with Nf = 4, 8 and 12 light or 
massless staggered fermions in the fundamental represen- 
tation. Fig. 1 collects our results for the scale-dependent 
mass anomalous dimensions 7 m (/i). For each of the three 
models, we calculate at least 1000 eigenmodes at several 
gauge couplings f3p on three to five lattice volumes as 
large as 32 3 x 64. 



Both of the top panels (a and b) of Fig. 1 refer to 
the Q CD-like 4- flavor system that we discuss in detail 
in Section III. In this system we are able to follow the 
evolution of the mass anomalous dimension from asymp- 
totic freedom (j m — >> 0) in the UV to the chirally broken 
IR. The Nf = 4 results for 7 m (/i) in panel (a) are plot- 
ted as functions of the lattice eigenvalue A. In Q CD-like 
systems, dimensionful quantities such as A implicitly de- 
pend on the lattice spacing a, which is determined by the 
gauge coupling. In panel (b) we show the same 4-flavor 
results with A expressed in terms of a uniform scale, a^A 
corresponding to f3p = 7.4. 

The uniform curve in panel (b) indicates that our 
Nf = 4 simulations are in the scaling regime of the gaus- 
sian fixed point at g 2 = 0. We compare our results for 
the energy-dependent mass anomalous dimension 7 m (/i) 
to one-loop perturbation theory. Our non-perturbative 
results break away from the (dashed) perturbative pre- 
diction around j m > 0.4, but still produce a universal 
curve well into the chirally broken regime, covering close 
to two orders of magnitude in energy scale. The results 
presented in Fig. lb show that the method to extract the 
anomalous dimension is well controlled and reliable even 
on very small physical volumes with vanishing fermion 
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masses. 

Panel (d) on the bottom right of Fig. 1 shows 7 m (/i) 
for the Nf — 12 system with various gauge couplings j3p. 
These results, discussed in Section IV, are very different 
from the Nf = 4 case. For all but the weakest coupling 
{Pf = 6.0), the anomalous dimension increases towards 
the UV, the opposite of the expected behavior near an 
asymptotically- free fixed point. In the infrared, all of our 
12-flavor results with different gauge couplings appear to 
approach a unique value as A decreases. We interpret this 
behavior as a strong indication of infrared conformality 
and identify the common A — » limit as the scheme- 
independent mass anomalous dimension 7^ = 0.32(3) 
that characterizes the conformal theory at the IR fixed 
point. 

Also unlike the 4- flavor case, the 12-flavor results for 
these gauge couplings cannot be rescaled to a unique 
curve. This is consistent with an irrelevant gauge cou- 
pling at the conformal IR fixed point, which does not 
allow a well-defined lattice scale [24]. It is also striking 
how much the 12-flavor j m changes with A at the stronger 
couplings Pf = 3.0 and 4.0. If one were to consider only 
a single gauge coupling and neglect the A — >• extrapola- 
tion, the resulting 7^ could be dramatically different [25] 

Finally, panel (c) on the bottom left shows the 8-flavor 
results that we also consider in Section IV. Again, these 
results do not indicate Q CD-like, asymptotically free UV 
behavior, but they also differ compared to Nf = 12. The 
anomalous dimension shows very little dependence on A, 
but changes with the gauge coupling. At our strongest 
accessible coupling ftp = 4.65, the anomalous dimension 
is 7m ^ 1, yet we do not observe spontaneous chiral 
symmetry breaking even on our largest volumes. We find 
Nf = 8 the hardest case to interpret, and we defer a 
detailed discussion of this system to a future publication 
that will also consider other observables, including finite- 
temperature transitions and the hadron spectrum. We 
conclude in Section V with some discussion of potential 
future improvements. 



II. MASS ANOMALOUS DIMENSION FROM 
THE MODE NUMBER 

A. Mode number scaling 

Conformal systems are chirally symmetric, so that the 
eigenvalue density p(X) vanishes at A = in the infinite 
volume, zero mass limit. The simplest scaling form valid 
for small A is p(X) oc A a , leading to the RG-invariant 
mode number 



'(A) = V f 



p(u)du oc VX 



l+a 



(2) 



Under a renormalization group transformation with scale 
factor s, the volume V — » s 4 V while A — > A/s 1+7 ™ in 
the infrared limit, where 7^ is the scheme-independent 
mass anomalous dimension at the IR fixed point. The 



renormalization group invariance of the mode number, 
VX 1+a = s 4 V (A/s 1+7 ^) 1+Q! , therefore relates the expo- 
nent a and 7 m as [11] 



1 



A -> 0. 



(3) 



The eigenvalues A define an energy scale: large A cor- 
respond to the UV while small A probe the infrared 
dynamics. By introducing an energy-dependent scal- 
ing exponent a(A), we can generalize the scaling form 
p(X) oc A a(A) to obtain 



l+7m(A) 



a(X) 



(4) 



This behavior is sketched in Fig. 2 for idealized (infinite- 
volume, zero-mass, continuum) IR-conformal and chi- 
rally broken systems. Considering only asymptotically 
free theories, j m — > as A — » 00 (the UV), which by 
Eq. 4 corresponds to a — ^ 3, reproducing the known scal- 
ing p(X) oc A 3 in free field theory [26]. For IR-conformal 
theories, 7 m — » 7 m as A — >• (the IR), to reproduce 
Eq. 3. In between these two extremes we obtain a scale- 
dependent exponent that connects the limiting UV and 
IR values, < 7 m (A) < 7 m . (We discuss the behavior of 
Fig. Id for N f = 12 in Section IV.) 

Chirally broken systems can be described similarly, as 
shown in the right panel of Fig. 2. The main difference is 
p(0) 7^ due to spontaneous chiral symmetry breaking, 
which does not follow the scaling form p(X) oc X a ^ x \ As a 
result, while naive application of Eq. 2 to chirally broken 
systems in the IR would produce a — >• and 7 m — » 3 as 
A — >• 0, this prediction has no physical significance. Such 
unphysically large 7 m (A) simply indicate the breakdown 
of the scaling form due to the onset of spontaneous chiral 
symmetry breaking (which is typically expected for j m > 

1.) 

In the context of lattice calculations, we are restricted 
to A smaller than the UV cutoff defined by the inverse 
lattice spacing a -1 (A < 1 in lattice units). While larger 
eigenvalues can easily be calculated, they are in a regime 
dominated by lattice artifacts where no universal behav- 
ior can be identified. 

To extract the scale-dependent exponents a (A) and 
7m (A) from the mode number, we simply perform a linear 
fit to the logarithms 

log [i/(A)] = (a(A) + 1) log [A] + constant, (5) 

using finite intervals in A and a jackknife analysis to de- 
termine uncertainties. All results we present here use 
0.015 < AA < 0.075, as listed in Tables I through III. 
For each ensemble, we choose AA by requiring that every 
linear fit considers at least 10 points. We find x 2 /dof < 1 
at most, and almost always much smaller. Changing A A 
does not significantly affect the central values of a(X). 
We keep the fit range as small as possible since we incor- 
porate AA into the error bands shown in Figs. 1 and 4, 
which are smeared out as AA increases. 
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FIG. 2. Cartoons of eigenvalue densities p(X) in IR-conformal (left) and chirally broken (right) continuum systems. In both 
cases, asymptotic freedom predicts j m —> as A — » oo in the UV. (On the lattice we are restricted to A smaller than the 
UV cutoff defined by the inverse lattice spacing a -1 .) In IR-conformal systems, j m —> 7m at the IR fixed point as A —> 0. 
Spontaneous chiral symmetry breaking produces p(0) oc (ipip) / 0, which does not follow the scaling form p(A) oc X a(yX \ 



B. Potential systematic effects 

The scaling form p(X) oc X a leading to Eq. 2 as- 
sumes that the system is in infinite volume with van- 
ishing fermion mass. Lattice calculations are necessarily 
carried out in a finite volume, and typically use non-zero 
fermion masses as well. Both finite volume and finite 
mass break conformal scale invariance, which can only be 
recovered by extrapolations to the infinite- volume, chiral 
limit. However, Ref. [21] found negligible finite- volume 
effects for the mode number measured on 24 3 x 64 and 
32 3 x 64 lattices, and also observed scaling behavior for 
surprisingly large fermion masses. 

While the results of Ref. [21] give us some confidence 
that these systematic effects are manageable, we must 
carry out our own tests to directly check the lattice re- 
sults presented in Sections III and IV. In particular, we 
use very small or vanishing sea fermion masses, and al- 
ways calculate eigenmodes of the massless Dirac oper- 
ator. We investigated finite-mass effects in Ref. [23], 
finding that light masses m < 0.005 produce indistin- 
guishable results for the spectral density p(X). While 
simulations of chirally broken systems require a non-zero 
fermion mass, the largest mass we use in this paper is 
m 0.0025. 

Because we use such small masses, the finite volume is 
a more serious issue. To address finite-volume effects, we 
carry out simulations with several different lattice vol- 
umes and gauge couplings, combining the results to ac- 
cess the infinite- volume physics. We start with a general 
review of the finite- volume effects on the spectral density 
p(A), extending our earlier investigations [23] to empha- 
size that these effects can be controlled even in the chiral 
limit and at weak gauge couplings. 

Fig. 3 shows p(X) for the Nf = 4 system on lattice 
volumes 24 3 x48, 16 3 x32 and 12 3 x24. In the left panel 
we consider the reasonably strong gauge coupling ftp = 
6.4, where the largest volume (with m = 0.0025) shows 



chiral symmetry breaking, p(0) ^ 0. (The ~30% drop in 
the smallest-A bin may suggest that the 24 3 x 48 volume 
is near the boundary of chiral restoration.) The other 
two systems are clearly volume- squeezed, and we observe 
a gap in the 12 3 x 24 eigenvalue density, which permits 
simulation in the m = chiral limit. While the small A 
region is affected by the finite lattice volume, this is only 
a transient effect. For A > 0.04 the two larger volumes 
are indistinguishable, and all three volumes converge to 
the same curve shortly thereafter. 

The right panel of Fig. 3 shows the Nf = 4 eigen- 
value density from the same lattice volumes, now with 
/3p = 7.0. This coupling is significantly weaker (the lat- 
tice spacing at /3p = 7.0 is approximately half of that 
at /3f = 6.4, Eq. 6), and we encounter no obstacle to 
working directly at m = 0: all three systems are volume- 
squeezed with gaps that grow smaller as the volume in- 
creases. Nevertheless at sufficiently large A > 0.15 the 
different p(X) again overlap, indicating infinite- volume 
results. As the two panels of Fig. 3 illustrate, it is pos- 
sible to identify infinite- volume behavior even when vol- 
umes in the p-regime are combined with strongly volume- 
squeezed systems, both with small and exactly vanishing 
fermion masses. 

The finite-volume effects we observe in p(X) at small A 
influence the fit (Eq. 5) that we use to extract the mass 
anomalous dimension. The left panel of Fig. 4 shows 
7m (A) for Nf — 4 with /3p = 6.4, which comes from the 
data in Fig. 3. On the 24 3 x 48 volume the anomalous 
dimension is large, 7 m (A) > 1, consistent with the chiral 
symmetry breaking established by Fig. 3. In volume- 
squeezed systems, the finite volume pushes j m — ^ as 
A — » 0. This behavior is also unphysical, and indicates 
the breakdown of the scaling form p(X) oc X a ^ due to 
finite- volume effects. Since we compare several lattice 
volumes in Fig. 4, we can easily identify these transient 
effects by observing where the results for a given volume 
break away from the universal curve. (We omit these 
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FIG. 3. Volume dependence of the Nf = 4 spectral density p(A), normalized per continuum flavor. We calculate 1000 
eigenmodes on each lattice volume 24 3 x48, 16 3 x32 and 12 3 x24 (Table I). Left: The stronger coupling ftp — 6.4 exhibits chiral 
symmetry breaking with p(0) > on 24 3 x48 lattices with m = 0.0025. Right: At the weaker coupling ftp — 7.0 we encounter 
no obstacle to working directly at m = on all three volumes. The insets enlarge the small-A behavior. 
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FIG. 4. Predictions for 7 m from the scaling of v with A (Eq. 2), illustrating finite- volume effects at small A. We measure 
1000 or 1500 eigenmodes on each lattice volume 24 3 x48, 16 3 x32 and 12 3 x24 (Tables I and II). Left: For Nf = 4 at coupling 
Pf — 6.4, we observe 7 m > 1 on 24 3 x48 lattices with m = 0.0025, consistent with chiral symmetry breaking (Fig. 3). Right: 
For Nf = 12 at @f = 4.0, we also include 18 3 x36 lattices with m = 0. 



transients in Fig. 1.) 



III. RESULTS FOR QCD-LIKE SYSTEMS 



Finally, the right panel of Fig. 4 considers the 12-flavor 
system with /3 F = 4.0 and m = 0.0025 on 24 3 x48, 16 3 x32 
and 12 3 x24 lattices, also including m = on 18 3 x36 lat- 
tices. Again, after the small-A transients, the 7 m from 
different volumes form a universal curve indicating that 
finite-volume effects are under control. In Fig. Id we see 
that at the weaker couplings ftp = 5.0 and 6.0, finite- 
volume effects increase, and even at large A the results 
from different volumes don't overlap perfectly. This sit- 
uation calls for infinite-volume extrapolations that are 
not yet feasible: controlled extrapolations require data 
for v(\) at larger A on the larger volumes. Directly mea- 
suring more eigenmodes is too computationally expensive 
to be practical, so we are developing code to stochasti- 
cally calculate the mode number, as in Refs. [5, 21]. We 
revisit this issue in Section V. 



In this section we describe our results for SU(3) lat- 
tice systems with Nf = 4 light or massless staggered 
fermions. This model exhibits QCD-like behavior, with 
spontaneous chiral symmetry breaking, confinement, and 
a running gauge coupling driven by the perturbative 
gaussian fixed point. Our 4-flavor tests verify the applica- 
bility of our method to predict the scale-dependent mass 
anomalous dimension, illustrate the benefits of combining 
different lattice volumes and gauge couplings, and con- 
firm the validity of results obtained from volumes much 
smaller than the confinement scale. 

In our lattice studies, we generate ensembles of gauge 
configurations using the hybrid Monte Carlo (HMC) al- 
gorithm with a second-order Omelyan integrator [27] ac- 
celerated by an additional heavy pseudofermion field [28] 
and multiple time scales [29]. We use nHYP-smeared 
staggered fermions [30, 31], with smearing parameters 
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TABLE I. Nf — 4 lattice ensembles used in Figs. 1, 3 and 4. 
For each ensemble specified by the volume, fermion mass, and 
gauge coupling /3f, we report the total molecular dynamics 
time units generated with the HMC algorithm, the number 
of configurations on which we measure at least 1000 (1500) 
eigenvalues, and the fit range AA. 



Volume 


Mass 


Pf 


Total MDTU 


# meas. 


AA 




0.0025 


6.4 


920 


32 (6) 


0.015 




0.0025 


6.6 


635 


26 


0.015 


24 3 x48 


0.0 


7.0 


800 


31 


0.0225 




0.0 


7.4 


790 


30 


0.0325 




0.0 


8.0 


1000 


40 (40) 


0.04 




0.0025 


6.4 


1365 


41 


0.0325 




0.0025 


6.6 


1125 


46 


0.0375 


16 3 x32 


0.0 


7.0 


750 


47 


0.05 




0.0 


7.4 


750 


47 


0.05 




0.0 


8.0 


1400 


41 


0.0525 




0.0 


6.4 


1000 


50 


0.055 




0.0 


6.6 


1000 


50 


0.07 


12 3 x24 


0.0 


7.0 


1000 


62 


0.07 




0.0 


7.4 


1000 


61 


0.07 




0.0 


8.0 


1840 


86 


0.075 



(0.5, 0.5, 0.4); the motivations for nHYP smearing in gen- 
eral and these smearing parameters in particular are 
reviewed in Ref. [10]. For the fermions we use peri- 
odic boundary conditions in spatial directions and anti- 
periodic boundary conditions in the temporal direction. 
Our gauge action includes an adjoint plaquette term with 
coefficient /3 a related to that of the fundamental pla- 
quette term by /3a = — 0.25/3 >. This negative adjoint 
plaquette term lets us avoid a well-known spurious ul- 
traviolet fixed point caused by lattice artifacts, and im- 
plies f3p = 12 / g 2 at the perturbative level. We measure 
1000 to 1500 eigenmodes on thermalized configurations 
separated by at least 10 molecular dynamics time units 
(MDTU). The lattice ensembles used in the 4- flavor anal- 
ysis are summarized in Table I. 

In addition to measuring eigenmodes on these ensem- 
bles, we use the Wilson flow to determine the lattice 
spacing as proposed by Refs. [32, 33]. At our strongest 
gauge coupling f3p = 6.4, where the lattice spacing is 
largest and finite- volume effects appear minimal, we find 
v^o/a = 1.42(1) and w /a = 1.94(2). While there is 
no precise determination of the scales y^o and wq for 
Nf = 4, from the 2+1-flavor results y^o = 0.145 fm and 
wq = 0.176 fm reported in Ref. [33] we estimate that our 
lattice spacing is around 0.1 fm at this coupling. This 
estimate is rough but allows our coupling values to be 
compared with those of other simulations that use differ- 
ent lattice actions. 

At our weaker couplings, finite- volume effects are more 
severe. Instead of measuring and wo/a on these 

ensembles, we find it is easier to calculate relative scales 
using a step-scaling approach like that described in 



Ref. [34]. By comparing the renormalized running cou- 
plings on 24 3 x48, 16 3 x32 and 12 3 x24 volumes, we predict 
the relative lattice spacings 

o<6a/cl7A = 2.84(3) a 6 . 6 /a 7A = 2.20(5) (6) 
&7.0/&7.4 = 1.45(3) &8.o/&7.4 = 0.60(4). 

Thus the physical size of our configurations changes 
about an order of magnitude between the 24 3 x 48 vol- 
ume at f3 F = 6.4 and the 12 3 x 24 volume at f3 F = 8.0. 
These are preliminary results and the errors in Eq. 6 are 
conservative, but they suffice for the analysis we consider 
here. 

In the previous section we discussed how different vol- 
umes can be combined to obtain volume-independent 
predictions for the mass anomalous dimension. This is 
illustrated in the left panel of Fig. 4 for the 4-flavor sys- 
tem with Pf = 6.4, and in Fig. la for all five gauge cou- 
plings we consider. To combine results for different /3p 
we rescale the lattice eigenvalues so that they are all 
expressed in terms of a uniform scale, ajA corresponding 
to (3f — 7.4. In the chirally symmetric regime we take 



ota 



l+7m(A/0 



(7) 



The scaling dimension l+7 m (A) appears because A scales 
in the same way as m. This is easy to understand by 
recalling that the massive Dirac operator is just D m = 
D + m. 

In the chirally broken regime the scaling form p(X) oc 
\ a W n o longer holds (Fig. 2), and results for 7 m de- 
termined from Eq. 4 are not physical. Therefore, when 
7 m > 1 we modify Eq. 7 to 



A, 



A, 



&7.4 

dp 



(8) 



Our choice of 7 m = 1 as the value at which we switch 
from Eq. 7 to Eq. 8 is motivated by the conventional 
wisdom that chiral symmetry breaking sets in at 7 m > 1 . 

Applying Eqs. 6-8 to the 4-flavor results in Fig. la pro- 
duces the universal curve shown in Fig. lb. (Since the 
uncertainties in Eq. 6 are preliminary, we do not incor- 
porate them into the error bands in Fig. lb, which are 
therefore underestimated.) Every volume and gauge cou- 
pling we consider can be combined to cover nearly two 
orders of magnitude in energy, from the chirally broken 
regime in the IR to the perturbative regime in the UV. 
The dashed line in Fig. lb is the one-loop perturbative 
prediction of the anomalous dimension, fitted with one 
parameter to match the scale at A7.4 = 0.8. 

The 4-flavor model provides robust tests of our pro- 
posal in the relatively familiar context of QCD-like sys- 
tems. We observe that the systematic effects discussed in 
the previous section are under control, justifying our use 
of volumes much smaller than the confinement scale. The 
universal curve we obtain after rescaling with Eqs. 6-8 
demonstrates the power of combining multiple volumes 
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and gauge couplings, and shows that finite-mass effects 
are negligible for m = 0.0025. These results increase our 
confidence in our findings for systems unlike QCD, to 
which we now turn. 



IV. RESULTS FOR SYSTEMS UNLIKE QCD 

In this section we consider two examples of sys- 
tems whose dynamics appear qualitatively different from 
QCD: the 12- and 8-flavor SU(3) models. In both cases 
we observe scaling that is not driven by the perturba- 
tive gaussian fixed point. Such behavior could be due 
to conformal infrared dynamics, but might also occur in 
chirally broken systems due to non-perturbative effects 
at strong gauge couplings. 



A. N f 12 

Several groups have investigated the 12-flavor lattice 
system using a variety of methods [10, 23, 35-58], arriv- 
ing at different conclusions regarding its infrared behav- 
ior (see Ref. [20] for a recent review). Our own investi- 
gations favor the existence of a conformal IR fixed point 
for Nf = 12, which we directly observe in Monte Carlo 
Renormalization Group (MCRG) studies [41, 46, 57]. 
In more comprehensive explorations of 12-flavor sys- 
tems with many different gauge couplings and fermion 
masses [10, 52], we observed an unusual lattice 
phase where the single site shift symmetry ("6' 4 ") of the 
lattice action is spontaneously broken (see also Ref. [53]). 
We restrict our present analysis of the Dirac eigenmodes 
to ftp > 3.0, weak enough to avoid the lattice phase. 
Our Nf = 12 simulations in this range of couplings, 
which include volumes as large as 32 3 x 64 and 40 3 x 20 
with masses m < 0.03, do not indicate spontaneous chiral 
symmetry breaking [10, 52]. 

Table II summarizes the 12-flavor lattice ensembles we 
use in this paper. As described in the previous sections, 
we consider multiple gauge couplings and combine results 
from several lattice volumes as large as 32 3 x 64. We use 
fermion mass m = 0.0025 on four of the five volumes, 
and m = on 18 3 x 36 lattices. In Ref. [23] we investi- 
gated Nf = 12 systems with larger masses m < 0.03, and 
verified the mass independence of p(X) for m < 0.005. 

Fig. Id presents our results for the mass anomalous 
dimension at four different gauge couplings. At each 
results from up to five different volumes overlap. In the 
right panel of Fig. 4 we zoom in on this overlap for f3p = 
4.0. The difference between the 4- and 12-flavor results 
is striking. While the Nf = 4 j m always decrease as 
the energy scale increases, as expected from asymptotic 
freedom, our Nf = 12 results at the stronger couplings 
Pf =3.0 and 4.0 show the opposite behavior, increasing 
towards the ultraviolet. At the weaker couplings ftp — 
5.0 and 6.0, the anomalous dimension is roughly invariant 
across an order-of-magnitude change in scale. 



TABLE II. Nf = 12 lattice ensembles used in Figs. 1, 4 and 
5, with columns as in Table I. 



\//~\ 1 1 1 V V 1 o 

V UlLLIIlc 


IvlcLob 


PF 


Total MTjTTT 






^9 3 V fU 
OZ A \)*-t 


u.uuzo 


o.u 


1 ^70 


1 1 

± ± 


01 ^ 
u.uio 




o 009^ 


*\ o 

o.u 


1 9^0 

1ZOU 


1 3 


01 ^ 

U.UIO 




o 009^; 
u.uuzo 


o 
o.u 


1 07^ 
1U i O 


/in 

4U 


O 01 ^ 
U.U10 


9/1 3 v/ jo 
z4 X 4o 


u.uuzo 


a n 

4.U 


1UUU 


OA fo r >\ 

Z4 yZZ J 


U.U10 




009^ 

u .uuzo 


O.U 


1 0DD 

-LUUU 


40 ^Q^l 

^±U ycf J 


01 ^ 

U.U-LO 




009^ 




U.U 


1 9^0 

1ZOU 


OU ^OU J 


099^ 
u.uzzo 




o o 

u.u 


o.u 


1 9^0 

1ZOU 


^9 
oz 


01 ^ 

U.UIO 


1 S 3 w 
lo X oO 


u.u 


a n 

4.U 


1ZUU 


ou 


n no 
u.uz 




o n 
u.u 


^ o 

o.u 


1 9^0 

1ZOU 


fi9 

UZ 


099^ 
u.uzzo 




o o 

u.u 



u.u 


1 9^0 

1ZOU 


^9 
oz 


09^ 
u.uzo 




o 009^; 
u.uuzo 


o 

O.U 


zuuu 


/in 

4U 


O 01 ^ 
U.U10 


1 fi3 OO 

1U X oz 


u.uuzo 


zt o 


you 


4tU yO j 


099^ 
U.UZZO 




0.0025 


5.0 


1020 


40 (6) 


0.03 (0.35) 




0.0025 


6.0 


1130 


24 (24) 


0.0325 




0.0025 


3.0 


2000 


40 


0.0225 


12 3 x24 


0.0025 


4.0 


550 


40 


0.0325 




0.0025 


5.0 


900 


40 


0.0425 




0.0025 


6.0 


850 


40 


0.045 



This scale dependence of 7 m is not consistent with in- 
frared dynamics driven only by the gaussian fixed point, 
and does not allow us to rescale results for different /3p to 
produce a combined universal curve. Undoubtedly as we 
approach the perturbative gaussian FP at very weak cou- 
pling, the 12-flavor system will exhibit asymptotic behav- 
ior like that we see for Nf = 4. We observe initial indica- 
tions of this at Pf = 6.0 in Fig. Id, and for /3p > 7.0 we 
have found that j m clearly decreases with increasing A. 
(While we have investigated 7.0 < fi F < 10.0 on 24 3 x48 
and smaller lattices, finite- volume effects become increas- 
ingly severe at such weak couplings; we would need lattice 
volumes of 32 3 x64 or larger to gain sufficient control over 
these effects for /3p > 7.0.) 

Our 12-flavor results in Fig. Id can be interpreted as 
IR-conformal behavior. For weaker couplings the gaus- 
sian fixed point drives the system towards the conformal 
IR fixed point, but around /3p = 5.0 the behavior of 7 m 
changes. At /3p = 3.0 and 4.0 the anomalous dimen- 
sion grows in the UV and it is tempting to state that 
these lattice systems are 'on the strong-coupling side of 
the IRFP, p F < However, recall that the IR fixed 
point is located in an infinite-dimensional space of lattice- 
action terms [24]. f3 F is defined through the projection 
of the IRFP onto the one-parameter space of /3p. While 
the existence of the conformal phase and its IR fixed 
point is a universal property, the location of the IRFP 
in the action-space, and therefore j3 Fl is scheme depen- 
dent. The lattice action provides a regularization, but 
still allows many different renormalization schemes to be 
explored. Each of these schemes can predict a different 
j3 Fl as we have observed in MCRG studies [41, 46, 57]. As 
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discussed in Section 6 of Ref. [59], one can define non- 
perturbative renormalization schemes through physical 
observables as well. We generally need two independent 
observables, whose dependence on the bare gauge cou- 
pling and fermion mass determine the bare renormaliza- 
tion group j3 and 7 functions. When the system is in the 
chiral limit, a single observable suffices. Our observable, 
the mode number in the chiral limit, suggests j3p ~ 5 as 
the projection of the IRFP. 

The IR-conformal interpretation of our 12-flavor re- 
sults demands that in the limit A — >• we obtain a uni- 
versal value 7^ for the scheme-independent anomalous 
dimension at the conformal IR fixed point. The finite 
volumes of our lattice systems prevent us from directly 
investigating A = 0. At the smallest A ~ (9(0.1) that we 
can access, the 7 m from different f3p vary over a wide 
range 0.2 < j m < 0.6. This is due to the slow running of 
the gauge coupling, as has been observed and discussed 
in similar models [14, 60]. Even though the dependence 
of 7 m on the gauge coupling dies out fairly slowly, the 
results for different f3p do approach a common value in 
the infrared, as expected for an IR fixed point at which 
the gauge coupling is irrelevant. 

By extrapolating A — > for each set of results with 
fixed /3p, we predict 7^ = 0.32(3). This value is con- 
sistent with 3- and 4-loop perturbative predictions in the 
MS scheme (7^ = 0.312 and 0.253, respectively [61]), and 
comparable to other recent lattice results for Nf = 12. 
Ref. [51] obtains 0.4 < 7^ < 0.5 from IR-conformal 
finite- volume scaling of spectral observables, considering 
two relatively weak gauge couplings and large fermion 
masses m > 0.04. Ref. [56] considers smaller masses 
0.006 < m < 0.035 at a single gauge coupling, but argues 
against the existence of an IR fixed point for Nf = 12 
on the grounds that finite- volume scaling of different ob- 
servables predicts different 0.2 < j m < 0.4. 

While the eigenvalues of the Dirac operator can re- 
veal a surprising amount of information, Fig. Id makes it 
clear that this approach has systematic effects that must 
be understood and addressed. Even in the 12-flavor sys- 
tem that we argue flows to a conformal fixed point in 
the infrared, the energy dependence of the anomalous di- 
mension can be significant. Analyses that do not check a 
range of energy scales risk obtaining apparently very pre- 
cise but actually incorrect results. Extrapolation to the 
infrared limit is necessary, and may significantly increase 
the numerical uncertainties. In addition, the slow run- 
ning of the coupling near the IR fixed point can make the 
results sensitive to /3p, even though the gauge coupling is 
irrelevant at the IRFP. Investigating several gauge cou- 
plings is important to control this systematic effect and 
confirm that universal results are obtained from extrap- 
olations to the IR limit. We overlooked both of these 
issues in our earlier 12-flavor study [10], where we re- 
ported a considerably larger value for 7^ based on data 
at a single (relatively strong) gauge coupling, ftp = 2.8. 



B. Nf — 8 

Two-loop perturbation theory predicts the existence 
of an IR fixed point for SU(3) gauge theories with Nf > 
8.05, implying that the 8-flavor model is just barely below 
the conformal window. Since this two-loop perturbative 
fixed point is at very strong coupling, higher-loop cor- 
rections could be significant [61, 62]. Analytic estimates 
based on the Dyson-Schwinger equation [63, 64] or a con- 
jectured thermal inequality [65] put the lower edge of the 

conformal window around N^ ^12, leaving the 8-flavor 
system well in the chirally broken regime. Lattice calcu- 
lations that study step scaling either with Schrodinger 
functional [35] or MCRG [38, 57] methods do not contra- 
dict this expectation. However, it is important to note 
that non-observation of an IRFP in a step scaling study is 
not sufficient to draw a conclusion. It is also necessary to 
show that the study probes strong enough couplings on 
large enough volumes that the system undergoes sponta- 
neous chiral symmetry breaking. Otherwise it is possible 
that investigations at stronger couplings would reveal an 
IR fixed point. For Nf = 8 this condition has not yet 
been satisfied [7, 10, 35, 37, 38, 52, 66-69]. 

Lattice simulations cannot necessarily reach arbitrar- 
ily strong couplings: lattice artifacts can induce first- 
order transitions, which separate strongly-coupled lattice 
phases from the weak-coupling phase where the contin- 
uum limit is defined. This issue is especially important 
for systems around the lower edge of the conformal win- 
dow, where very strong couplings may be required to dis- 
tinguish IR conformality from chirally-broken dynamics. 
The Schrodinger functional calculation of Ref. [35] using 
unimproved staggered fermions encountered a clear first- 
order transition. Improving the action by smearing the 
fermions allows us to reach stronger couplings before ob- 
serving a different first-order transition that borders the 
JS^ lattice phase [10, 52]. As for Nf = 12, in this paper we 
consider only couplings weak enough to avoid the 8-flavor 
phase, though our strongest /3p = 4.65 is close to this 
transition. Table III summarizes these ensembles. While 
we have not been able to establish spontaneous chiral 
symmetry breaking in this range of couplings, even on 
volumes as large as 32 3 x 64 and 40 3 x 20 with masses 
m < 0.02 [52], neither did we observe an IR fixed point 
with the MCRG method [38, 57]. We continue to inves- 
tigate this puzzle and will report further results in future 
publications. 

We show our results for the 8-flavor anomalous dimen- 
sion in panel (c) of Fig. 1. They differ from both the 4- 
and 12-flavor cases. At each fixed coupling, we find 7 m to 
be roughly constant over the order-of-magnitude change 
in scale accessible from combining lattice volumes, with 
a slight tendency to increase both towards the UV and 
towards the IR. 7 m (A) increasing with larger A suggests 
that these systems are not in the basin of attraction of 
the perturbative fixed point. (At j3p = 7.0 we observe 
the system approaching the gaussian FP; while we have 
investigated 6.0 < /3p < 8.0 on 24 3 x 48 and smaller 
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TABLE III. N f = 8 lattice ensembles used in Fig. 1, with 
columns as in Table I. 



Volume 


Mass 


0F 


Total MDTU 


^ meas. 


AA 




0.0 


4.65 


224 


11 


0.015 




0.0 


4.7 


385 


24 


0.015 


24 3 x 48 


0.0 


4.8 


540 


37 


0.015 




0.0 


5.0 


435 


34 


0.015 




0.0 


5.4 


690 


52 


0.015 




0.0 


4.8 


960 


50 


0.015 


18 3 x36 


0.0 


5.0 


930 


52 


0.02 




0.0 


5.4 


1000 


50 


0.0225 




0.0 


4.65 


980 


25 


0.0175 




0.0 


4.7 


1250 


70 


0.0175 


16 3 x32 


0.0 


4.8 


595 


29 


0.02 




0.0 


5.0 


690 


39 


0.0225 




0.0 


5.4 


940 


36 


0.0275 




0.0 


4.65 


750 


47 


0.0275 




0.0 


4.7 


1250 


67 


0.0325 


12 3 x24 


0.0 


4.8 


1250 


87 


0.035 




0.0 


5.0 


1250 


87 


0.035 




0.0 


5.4 


1250 


43 


0.045 



lattices, we would need larger volumes to control finite- 
volume effects for f3p > 6.0.) The non-monotonicity of 
7m (A) prevents us from rescaling \p to obtain a combined 
universal curve. The A dependence of the 8-flavor anoma- 
lous dimension most closely resembles that observed for 
Nf = 12 at /3p = 5.0. As we move to stronger cou- 
plings, 7 m increases significantly but the overall structure 
of 7 m (A) does not change. 

These results, like our other investigations of the 8- 
flavor system, are hard to interpret. At a minimum, we 
observe that the anomalous dimension is large, 7 m ~ 1 
across a wide range of scales. We also see that differ- 
ent gauge couplings produce greater changes in j m than 
does evolution over an order of magnitude in energy. This 
behavior is what one would expect from approximately- 
conformal "walking" dynamics, but is also consistent 
with a strongly-coupled IR fixed point where the gauge 
coupling is irrelevant but runs slowly. 

Considerably more work is needed to reach more con- 
crete conclusions about the IR dynamics of the 8-flavor 
model. Our final observation here regards the JS^ phase 
that prevents us from investigating stronger couplings. 
We mentioned above that j3p = 4.65 is close to the tran- 
sition into this lattice phase, and we note from Fig. lc 
that at this coupling the anomalous dimension is 7 m « 1 
for all accessible A. Similarly, for Nf = 12 we encounter 
the JO* phase at (3p « 2.65, where 7 m ~ 1 for A > 0.1 
(consider the results for /3p = 3.0 in Fig. Id). This may 
be a coincidence, but we intend to explore whether there 
is a relation between the lattice phase and a large 
anomalous dimension in the ultraviolet. Such a connec- 
tion could be consistent with the interpretation that the 



$fi phase arises from a new relevant operator, such as 
the next-to-nearest neighbor interactions considered by 
Ref. [53]. A new relevant interaction implies the exis- 
tence of a UV fixed point where a continuum limit can 
be defined. Since we observe phase separation between 
the phase and the weak-coupling phase, it is likely 
that such a continuum limit would be non-trivial. We 
may speculate that such a UV fixed point might be asso- 
ciated with conformal infrared dynamics [62, 70], which 
merits further investigation. 



V. CONCLUSION 

We have shown how to extract the scale-dependent 
mass anomalous dimension 7 m (/i) from the renormaliza- 
tion group invariant Dirac operator mode number ^(A). 
We tested our method with 4-, 8- and 12-flavor SU(3) 
lattice gauge theories, to confirm that we control all sys- 
tematic effects and to demonstrate that by considering 
multiple lattice volumes and gauge couplings we can de- 
termine the anomalous dimension across a wide range of 
scales. 

In our numerical calculations we used nHYP-smeared 
staggered fermions and generated gauge configurations 
at very small or vanishing fermion masses. We measured 
1000 to 1500 eigenmodes of the massless Dirac opera- 
tor on several lattice volumes as large as 32 3 x 64. By 
combining different lattice volumes at fixed gauge cou- 
pling we can identify finite- volume effects and determine 
infinite- volume results in an energy range covering about 
an order of magnitude. 

We also consider many gauge couplings, which in the 
case of the 4-flavor system we are able to combine by 
rescaling with the lattice spacing. We predict a universal 
curve for the Nf = 4 anomalous dimension, which is con- 
sistent with one-loop perturbation theory in the ultravio- 
let. In the infrared we observe chiral symmetry breaking 
as expected for a QCD-like system. Our 4-flavor results 
thus demonstrate the validity of our method. 

Our 12-flavor results are very different from the 4-flavor 
case, and indicate the existence of an infrared fixed point. 
At strong couplings we observe j m increasing towards the 
ultraviolet, indicating that the system is not in the basin 
of attraction of the gaussian fixed point. Our results at 
different j3p cannot be combined to predict a universal 
curve. This implies that the concept of the ^-dependent 
lattice spacing is not applicable in the 12-flavor system, 
consistent with an IRFP where the gauge coupling is ir- 
relevant. Even though f3p is irrelevant, it runs slowly, and 
we observe that the dependence of our 12-flavor results on 
the coupling dies out fairly slowly. The universal anoma- 
lous dimension 7^ of the conformal fixed point can only 
be found by extrapolating to the A = infrared limit. 
We predict 7^ = 0.32(3), where the error is dominated 
by the uncertainties of A — » extrapolations at different 
3 < p F < 6. 

Our results for the 8-flavor model also illustrate the ef- 
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FIG. 5. Finite- volume effects on the mass anomalous dimen- 
sion predicted from the scaling of v with A (Eq. 2). Results 
are for Nf = 12 as in Fig. 4, except at the weaker coupling 
0f = 5.O. 



fects of the slowly-running gauge coupling. We find that 
7 m does not change significantly over the range of scales 
accessible at fixed /3p, but shows a strong dependence 
on the coupling itself. This behavior is consistent with 
walking dynamics, but it is also possible that the Nf = 8 
system is IR conformal with a fixed point at very strong 
coupling. More work is required to either observe this 
IRFP or rule it out by demonstrating spontaneous chiral 
symmetry breaking. Regardless of the ultimate IR fate 
of the 8-flavor model, our observation that the anoma- 
lous dimension is large over a wide range of energy scales 
makes this a very interesting system to study. 

The main limitation of our present work is the number 
of eigenvalues we can calculate on the larger lattice vol- 
umes. For example, in Fig. 1 there are some signs that 
for our weaker gauge couplings the results obtained on 
different volumes don't perfectly overlap. To provide a 
clearer illustration of this issue, in Fig. 5 we zoom in on 
the 12-flavor 7 m at ftp = 5.0. This figure suggests that at 



fixed A the predicted 7 m for this gauge coupling system- 
atically decrease as the volume increases. However, even 
measuring 1500 eigenmodes does not allow us to perform 
controlled extrapolations to the infinite- volume limit. 

Stochastic estimation of the mode number, as proposed 
in Ref. [5] and used by Ref. [21], is a promising way to 
increase the reach of our approach. We are currently 
developing code to perform this calculation. However, 
this approach is made more difficult by the requirement 
that we investigate a range of energy scales up to A < 0.5. 
In each fit range with 0.015 < AA < 0.075, we need to 
calculate v(X) for ~10 values of A. This is an expensive 
proposition (consider that Ref. [21] only carried out a 
few dozen such measurements in total), but is still more 
efficient than direct calculations of more eigenmodes. We 
look forward to discovering how the stochastic approach 
will perform, and what physics results it will provide. 
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